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Abstract 



A recent analysis of the 4-point correlation function of the passive scalar advected by a time- 
decorrelated random flow is extended to the A^-point case. It is shown that all stationary- 
state inertial-range correlations are dominated by homogeneous zero modes of singular op- 
erators describing their evolution. We compute analytically the zero modes governing the 
A^-point structure functions and the anomalous dimensions corresponding to them to the 
linear order in the scaling exponent of the 2-point function of the advecting velocity field. 
The implications of these calculations for the dissipation correlations are discussed. 
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1 Introduction 



There has been much effort lately to understand the behavior of a scalar quantity passively 
advected by a random flow with a Gaussian statistics decorrelated in time This simple 
model, of its own interest, has served as a prototype of a turbulent system. It is believed 
that its behavior may teach us important lessons about the fully developed hydrodynamical 
turbulence. One of the interesting aspects of the passive scalar which has been recently 
understood p|, ^, ^ is the origin of the breakdown of Kolmogorov inertial-range scaling in the 
higher structure functions of the scalar. It has been realized that the dominant contribution 
to the structure functions comes from the zero modes of the differential operators describing 
the stochastic evolution of the correlation functions of the scalar. In this note, we extend 
the results of ref. by presenting the computation of the anomalous dimensions of the Ap- 
point structure functions in the first order of the parameter ^. Exponent ^, which in |0 was 
denoted k and in 2 — 7, is the growth rate of the 2-point structure function of the velocities 
of the advecting flow. The present work was motivated by where a similar analysis in the 
first order in inverse dimension was sketched. 

The equation governing the passive scalar in a turbulent flow is: 

dtT + {u-V)T-iyAT = f. (1) 

Here T{x,t) describes the scalar, e.g. the temperature, and / the forcing term whose role is 
to compensate the dissipation caused by the term proportional to the molecular diffusivity u. 
The velocity field u with V • m = is supposed to be random. We shall work in d > 3 space 
dimensions and shall assume homogeneity, isotropy and parity invariance of the advecting 
flow and of the forcing. 

The statistics of the forcing term is assumed to be Gaussian with mean zero and 2-point 
function 

{f(x,t)fiy,t')) = CC-j^)5it-t') (2) 

The rotation- invariant function C{x/L), which could be chosen to be a Gaussian, varies on 
scale L. 

The statistics of the velocity field, independent of the forcing, is also supposed to be 
Gaussian with zero mean and with the 2-point functions 

(M"(a;, t)u'^{y, t')) = D^'^^x - y)6{t - t') with d^D"^ = 0. (3) 

To analyze the scaling property of the scalar correlation functions we shall use the following 
expression for D"^{x) = D(0)5"^ - with 

r^{x) = D (^(rf + e - 1) - \xf (4) 

where is a parameter, < ^ < 2. Clearly, this distribution for u is far from realistic. It 
mimics however the growth of the correlations of velocity differences with separation distance, 
typical for turbulent flows. The fact that the 2-point functions (0,^ are white noise in time 
is crucial for the solvability of the model. The parameter ^ fixes the naive dimensions under 
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the rescalings x — *• fix, L — »• fiL. The naive dimension of u is ^/2 and of T is (2 — ^)/2. 

( \^/? 

Scale L serves as an infrared cutoff and the "Kolmogorov scale" "<] = \^] as an ultraviolet 
cutoff. 

We shall be interested in the correlation functions of the scalar in the inertial range 
?7 -C x -C L. The main result of this note is that in this range the stationary-state, equal- 
time, even structure functions scale with the anomalous exponents pat as 

((T(a;,t)-T(0,t)n = aJ^^ \xf'^^''l^ (5) 



\x\ 



with 

The exponents are universal depending only on ^ but the amplitudes a^q are not: they depend 
on the shape of the covariance C . The error term is bounded by (9((L/|a;|)~^"'"'^*^^^|a:|*^^~^^^/^) 
so it is strongly suppressed for large As it should be, the pat's satisfy the Holder 

inequality > ^^p^. More precise descriptions and statements will be given below. The 
formula (H) agrees with the N = 4 result of and with the —expansion of [§, If. 

Following ref. P], we shall derive the values of the anomalous exponents by analyzing in 
perturbation theory in ^ the zero modes of differential operators characterizing the station- 
ary state. Although for ^ = one observes a purely diffusive behavior of T and for ^ > 
an inertial energy cascade, the zero modes differ little in both cases, the different physics 
arising from their cumulative effect. As already stressed in [0, this resembles the situation 
in the renormalization group analysis in field theory or statistical mechanics where relevant 
perturbations, controlable in the single scale problem, may have large effects on the behavior 
of the system. As in the renormalization group study of critical models, the 1*^* order pertur- 
bative corrections to the zero modes lead to the resummation of leading infrared logarithms 
in the perturbation expansion of the structure functions in powers of ^. Pursuing further 
the analogy we shall introduce, as in the perturbative renormalization group, the notion of 
matrix of anomalous dimensions, see also 0. 

Our results can be used to deduce the scaling properties of the correlation functions of 
the dissipation field, which we denote by e(x), as discussed for example in 0, |[. At finite 
diffusivity z/ 7^ 0, the dissipation field is defined by e(x) = ulimx'-^x (VT)(a;') ■ (VT){x). This 
is a sensible definition since at finite u the correlations of T and their first derivatives are 
not singular at coinciding points (the higher derivatives are). In the limit u ^ 0, we have 
alternative definitions: 

e{x) = lim u lim (VT)(a;') • (VT)(a;) (7) 

>0 x'^X 

or 

e(x) = lim Ud''^(x-x')dxc.d^,A lim T(x')T(x). (8) 

x'^x ^ ^ ^ i/^O 

The order of the limits in the first definition is crucial since when u ^ and for small |x — x'l, 
T{x)T{x') ~ I a; — modulo more regular terms so that (VT) (a;) ■ (VT) (a;') ~ \x — x'\~^ 
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and becomes singular. The non-commutativity of the hmits p ^ and x' ^ x is at the 
origin of the dissipative anomaly. The second definition of e(x) is in the spirit of the operator 
product expansion in the i/ = theory. Using the Hopf identities (p!OD for the correlation 
functions, we shall argue that both expressions for the dissipation field e{x) coincide for 
^ < 1. The mean dissipation e = (e(x)) is equal to |C(0) i.e. to the mean injection rate of 
energy. The dissipation field has zero naive scaling dimension since T^(x) and d{x)V'^ have 
opposite naive dimensions. However, as a consequence of the relation (^, one finds that e{x) 
acquires an anomalous scaling. In fact, the definition (|) and eq. ([T6| ) allow to compute any 
structure functions with (non- coincident) insertions of the dissipation field. For example, 
the connected 2-point function of e scales as 

(<^),<0))^~ (^^y (9) 

and it decreases with in agreement with the physical picture of the dissipation being 
a local process. Similarly, the n-point functions of e scale with exponents p2n- The short 
distance singularity in eq. is an unphysical artifact of the assumed short distance scaling 
of the advecting velocity, mollified in real systems by viscosity. 

The same method allows to obtain information about the dissipative terms appearing in 
the differential equations obeyed by the structure functions and to compare our results with 
the early attempts to calculate the anomalous exponents of the passive scalar and with 
the more recent ideas ^ about the behavior of the probability distribution functions in the 
turbulent systems. 



2 Inertial range scaling and the zero modes 

The correlation functions of T satisfy the (Hopf) identities which may be deduced using 
standard functional manipulations of stochastic differential equations, see e.g. or, for the 
present context, |]TD[. In the stationary state, the odd correlations vanish and the even ones 
satisfy at equal times the identities 

J2 -^A, + '-D{d- 1) Mn (T(xi) . . . T{xn)) 
i=i ^ ^ 

= j:C{x,k/L) (T(xi).^..,.r(x^)) (10) 

j<k 3 k 

with Xjk = Xj — Xk , Aj denoting the Laplacian in the Xj variable and with M. n standing for 
the differential operators given by 

N 

\D{d-l)MM = -^{Y.'^.,? + \Y. d-^{x,,)d^.d^,^. (11) 

The first operator on the r.h.s. of eq. (|T^) is zero by translation invariance and A^at is a sum 
of the 2-body operators. For u > 0, the operators appearing on the l.h.s. of equations ( p!OD are 
elliptic and positive. We may use their Green functions to solve the equations inductively. 
This will produce equal-time stationary correlators decaying at infinity. Physically, they 
describe the stationary state obtained by starting e.g. from a fixed localized configuration 
of the scalar and waiting long enough. 
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Notice that at ^ = the operator Ai n reduces in the translation invariant sector to the 
Laplacian in N variables Xj: A4n\^=o = — A^r = — Z^^i ^j- This implies that T becomes a 
Gaussian field at ^ = with the higher correlation functions built in the standard way from 
the 2-point ones. The stationary state coincides then with that of the forced diffusion with 
the effective diffusion constant equal to z/ + ^D{d — 1). 

We shall describe the inertial range correlators by taking the limit u at fixed positions 
Xj and fixed large infrared cutoff L. It is not important that positions Xj be disjoint as long 
as we do not take derivatives of the correlators, see the remarks after eq. (^. In the limit 
1/^0, the correlation functions satisfy eq. (p!0[) but without the terms z/A. These equations 
completely determine the inertial range correlators up to zero modes of operators 
Physically, the zero mode contributions are fixed by the fact that we consider the system 
which is the limit of the one with positive diffusivity u. Mathematically, this means that 
in order to inductively solve eqs. (|10|) we should use Green functions of the singular elliptic 
operators Such Green functions are limits of the Green functions of the non-singular 

operators corresponding to the u > case. It has been argued in ref. 0, |^, Q that the 
zero modes of operators Ai n effectively appear in the inertial range correlators and give the 
dominant contributions in the limit L oo. 

The zero modes in question are homogeneous under dilation, invariant by translations, 
rotations, parity and symmetric under permutations of points. Since at is a sum of 
2-body differential operators, zero modes of Ai n-i lead by symmetrization to zero modes of 
Mn. More precisely, if /Ar_i(a;i, ■ ■ ■ ,X7v-i) is a zero mode of TWtv-i, then 

fN{xi,...,XN) = XI • • ■ '^<^{A^-1))' (12) 

(Tg5jv 

where the sum is over the permutations of objects, is a symmetric zero mode of Ai^- 
These zero modes will never contribute to the structure functions {Y[j (T{xj) — T{yj))). At 
= 0, the zero modes of A^at = — A^r are polynomials. For any even N > 2 there is only 
one "new" zero mode of scaling dimension A^ that cannot be expressed as a symmetrized 
sum of the zero modes of J^n-i- We shall denote it by Eq (of course, Eq is defined only up 
to a combination of the latter). Explicitly, 

Eo{x„...,x^)= Yl n ^i-i^ + [•••] (13) 

pairings 

l<i_ <i+<iV 

where the dots [ . . . ] refer to quantities which may be written as a (symmetrized) sum of 
functions depending only on A^ — 1 variables. 

The 2-point function of the scalar in the inertial range is ^ |10| 

(r(xi)T(x2)) = const. - (2_^)^'^(^_i) l^i2p"^ + OiL-'\x,,\'-^) (14) 

with const. = 0{L'^^^). It follows that at ^ = 0, where T becomes a Gaussian field, 

(T(xi)...T(xjv))k=o = c%Eo{xi,...,xn) + [..•], (15) 

where = ^ ■ "^^^ error not contained in the [ . . . ] terms is bounded by 

0(L-2(max|x,-fc 1)^+2). 
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Upon switching on positive ^, the symmetric zero modes of degree N will evolve to 
zero modes of A^jv with a modified homogeneity. They may be found by the degenerate 
perturbation expansion. Again, only one of them will not come from the zero modes of 
J^N-i- We shall call it Fq. Although for ^ positive, T is no longer a Gaussian field, its 
correlation functions may be inductively computed from eq. (|10|). In particular it is easy to 
see that the simple expressions 

^<j<k<N 

where the coefficients An = ^(7^/2)!''' ( (2-g)"D'(d-i) ) ^ H ('^+(2 — 00^^5 satisfy the version 

of eq. (0) with u = and L = 00. This scaling solution obviously leads to vanishing higher 
structure functions and cannot give the right answer for the inertial range correlators. The 
homogeneous zero modes of the operators Mn, which enter already at the first inductive 
step (the constant in eq. (Q)), modify the answer. At further inductive steps, the previous 
step modifications will induce new ones which, however, all give rise to combinations [ . . . ] 
of functions depending on fewer variables except, eventually, for the terms proportional to 
zero modes of Ai^. If the homogeneity degree of the zero mode is smaller than (2 — ^)y, 
the proportionality constant will contain a compensating positive power of L and may give 
the contribution dominating the large L structure functions if the zero mode is not of the 
[ . . . ] type. Indeed, for small positive ^ there is only one non-[ . . . ] zero mode which we 
have denoted by Fq. Its homogeneity degree is {2 — — Pn with positive p^, as will be 
demonstrated below. 

As a result, for small > 0, 

(T(xi)...r(x^)) = cnL''^Fo{xi,...,xn) + [...] (16) 

with the non-[ . . . ] error bounded by 0{L^'^^'-^^^\max\xjk\)^'^'^^^^^^)- For ^ not very small, 
the perturbations of zero modes which at ^ = have degree higher than N may eventually 
enter the interval of scaling dimensions smaller than (2 — and give non- negligible or 
even dominant contributions to the structure functions. The large L and ^ — limits of the 
correlation functions of T do not commute since the terms scaling with different powers of 
L become degenerate for ^ = 0, see 0. These limits, however, do commute for the structure 
functions involving only the Fq contribution scaling as L^*^^-' and the error bounded by 
Fq\^=q = Eq, it follows by comparison of ([T5|) and (|16D that the amplitude 
cn = c% + 0{^). The 0(0 contributions to the amplitudes cn depend on the shape of 
covariance C and hence are not universal. 

The relation (|T6]) implies the behavior (^ of the A^-point structure functions Sn{x) = 
{ {T{x) - T(0))^). In the Gaussian limit, 

SN{x)k=o = (17) 

where = (^y/^ \ Dd{d-i) ) ■ follows from the continuity of the structure functions at 
^ = that the amplitude a at in eq. (H) is equal to + 0(0. 

In the perturbation expansion in powers of ^, 

FQ = EQ + ^GQ + 0{e). (18) 
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In the next sections, we shall compute the 0{^) contribution Go (modulo [ . . . ] terms). 
Inserting the decomposition ([ISD into ([T6|), we obtain an asymptotic expression for the struc- 
ture functions which, although obtained by the first order zero-mode analysis, contains all 
orders in ^ resumming the series J^'^n^^i^ogL)"' of logarithmic infrared divergences appear- 
ing in the expansion of the structure functions in powers of ^. This is the situation well 
known from the perturbative renormalization group where the first order approximation to 
the single renormalization group step leads upon iterations to the resummation of the leading 
logarithms in the perturbative expansion of correlation functions. 



3 Anomalous dimensions at 0{^) 

Let us discuss the perturbative calculation of the homogeneous zero modes of A4n- At the 
first order in ^ we have 

MN = -AN + ^VN + 0{e) 
with Atv the Laplacian in A^- variables and Vn given by 

Vn= i^"' log l^.^l - 71^^f^] ^^x" - 77^^^ • (19) 

Note that, since A^at is a homogeneous operator of dimension ^ — 2, we have 

[E % ' Vn] = - - 2Vn . (20) 



Let E' be a symmetric homogeneous zero mode of AiN\^=o of degree A^. We shall search 
for the zero mode of 7W at of the form F = E + + 0(^^). The zero-mode equation gives 
at the order linear in ^ 

-AnG + VnE = 0. (21) 

The solutions G of this equation are clearly defined up to zero modes of A^v. Note that due 
to the scaling properties of E and V^, 

- An (E^^i^ -n)g = - -N + 2)AnG 

= -{j2 % - + 2) VnE = AnE = . (22) 

Hence the function E' = (j2x'jdx" — G is necessarily a zero mode of A^r. We shall show 
that there exist solutions G of eq. (^) such that E' are homogeneous polynomials of degree 
N. Such solutions are defined up to degree N zero modes of Atv but this ambiguity does 
not show up in E'. We obtain this way a linear transformation 

T -.E^ E' 



of the space of symmetric homogeneous zero modes of A at of degree N. If (Ea) is a basis of 
this space then the matrix (F"^) of this transformation given by E^ = T^-fjEa plays the role of 
the matrix of anomalous dimensions at first order in ^. Indeed, if E = v^Ei, is an eigenvalue 
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A eigenvector of the transformation F, i.e. if {v^) is an eigenvector of matrix (F"^), then, for 
the corresponding solution of eq. (|2l|) , we obtain 



(23) 



or 



(E -N-^X){E + ^G) = 0{e) 



(24) 



which means that E + C,G is homogeneous of order + up to 0(^^). Hence the ho- 
mogeneous zero modes of A^at are perturbations of the ^ = zero modes corresponding 
to eigenvectors of the matrix of anomalous dimensions. If the matrix (F°^) is not totally 
diagonalizable then there will be logarithmic corrections to the zero- mode homogeneity ||10|| . 

Reflecting the fact that all but one zero modes of A^a? are obtainable from those of A4 n-i 
by symmetrization, the matrix of anomalous dimension is block triangular. Namely: 



if the matrix is written in a basis {Eq, {Ea')) where Eq is the zero mode defined in eq. ([I3D 
and {Ea') forms a basis of the degree zero modes arising by symmetrization of functions 
depending on at most — 1 variables. The matrix element F°q is necessarily an eigenvalue of 
the matrix (F*^^) . If by an adequate choice of the [ . . . ] terms in its definition Eo becomes the 
corresponding eigenvector of the transformation F then Fq = Eq + C,Go + 0{C,^) describes the 
perturbed homogeneous zero mode of A^tv and F°q gives the C(0 correction to the scaling 
exponent of Fq equal to N in the leading order. The anomalous exponent of the A^-point 
structure function is therefore given by 



In principle it may happen (although it does not at least for A^ = 2, 4, 6) that F'^q is a 
degenerate eigenvalue of F and there is no corresponding eigenvector Eq. It is easy to see, 
however, that even in this case there exists a zero mode Fq = Eq + ^Gq + 0{^^) of A^tv which 
is homogeneous of degree A^ + ^ F°q -|- 0(^^) up to [ ... ] terms (the homogeneous terms 
are accompanied by the ones with powers of logarithms, the latter appearing in the [ • • • ] 
subspace). Such modifications would not effect the analysis of the structure functions. 




PN = -e(y + n) 



(25) 



since the naive A^-point scaling dimension is {2 — ^)y- 



4 Leading order corrections to the zero modes 



Let us return to the analysis of eq. (^Tj). We shall search for the solution G in the form 




(26) 



with Hjk and H polynomials of degree A^. For such a solution. 




(27) 
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would necessarily be a zero mode of A at of degree A^, as required. We shall see that there 
indeed exist solutions of ( ^1]) of the form (|^) (unique up to degree N zero modes of A^r) and 
that the polynomials Hji^ scale as \xjk\^ when \xjk\ —>■ assuring that the logarithms in ( |26D 



do not lead to divergent singularities in the correlation functions of T at coinciding points. 
Note, however, that the divergences at coinciding points start to appear in the correlation 
functions involving double derivatives of T or products of two first derivatives. 

The substitution of the Ansatz (^) into eq. (^) gives a set of three equations for Hj^ 
and H: 

^NHjk = Vj-VkE, (28) 
(^d-2 + x,,-V,k)H,, + ^^^^{xJ,x%d,.d^,)E = -\x%K,,, (29) 



^nH = Y.K,k (30) 



where Kjk are polynomials of degree N — 2 and Vjfc = Vj — Vfc with Vj = {dx'^)- Eq. (|30D 
is for free since any polynomial of degree — 2 is in the image of A^r acting on polynomials 
of degree A^. Thus, given the solution Kj^ of eqs. (^,|^), there always exists a degree A^ 
polynomial H solving eq. (PU[) and it is unique up to the zero modes of Ajy. 



We are thus left with solving eqs. (p8| , p9D . We shall first prove that there is a unique 
solution of these equations and then we shall produce the solution when the initial zero 
mode E is the "new" zero mode Eq defined by eq. (0). Notice that this is now eq. ( pSj) 
which implies that Y^j^kHjk is a zero mode of Aat. By symmetry, we may specialize eqs. 
(pSj p9D to j = 1, k = 2. Let us work in the variables x = ^^^^^ , y = Xi2 and X3, . . . , xjy. We 



have: X12 ■ V12 = 2?/ ■ V^^. Since y ■ Vy counts the degree in y, we shall decompose all terms 

Trip) 

^p=o -"12 ) 



entering in eq. (p9D into a sum of terms of given degree in y. Namely: H12 = J2p=o 



Ku = E^To' and 



2{d 



1 ^ 



p=2 



with h[1\ and E^p) homo geneous polynomials in y of degree p. The operator {d — 2 + 
2y ■ Vy) is invertible on such homogeneous polynomials. Eq. (^) implies then that 

= ^3Y^(^^"^-y^ir'0 2<p<N. 

The fact that i/j^'' = -^^12^ = implies that Hu scales like {xul"^ when xi — > X2- Equation 
may then be rewritten as 



2AyHil^ = (Vi ■ V2^)(^-') - A^ifJr'^ 



N 



for p > 2, where A^ = ^-A^. + Aj. With the use of the previous relation between 
and the latter equation takes the form 



Ay(y'K 
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for some recursively known homogeneous polynomials fp. These equations may be solved for 
since A^y^ is an invertible operator on the space of homogeneous polynomials of fixed 
degree. 



Let us find the deformed zero mode Fq which at ^ = reduces to Eq of eq. (|T3|). To solve 
eqs. (^,^9D (by symmetry, we may again set j = 1 and k = 2), we first have to compute 
(Vi ■ V2)E and [x'^^x^^2dx'id^^)E. 

(Vi-V2)i? = -2id + N-2)Y:' n + (31) 

{x'^,x^,,d..d^,)E = -2xi,y! n ^li^ 

(i-U) 

+ 2 y: i4,-4jK4k~4k)Y." n ^i-i^ + (32) 

3<j<k<N 

where J2' denotes the sum over pairings {(/_,/+)} with 3 < Z_ < /+ < and J2" a similar 
sum but with, additionally, l± ^ j, k. The symbol [ . . . ]i2 refers to a sum of terms which do 
not depend on at least one Xp with p > 3. Recall from the proof of existence of solutions of 
eqs. ( p8| , p9D that H12 has to scale at least as |a;i2p as Xi X2- It follows then that it must 
be of the form 

ifi2 = aLu + b xl,J2' n ^1-1+ + [ • • • ]i2 (33) 

for some coefficients a and b where 

Lu = E i^l - - 4,) E" n ^i-i, ■ (34) 

3<j<k<N 



Note two properties of L 



12 • 

(xia- V12) L12 = 4Li2, (35) 

AL12 = -4(iV-2)E' n + [■■■h- (36) 

(i-U) 

The first relation just means that Lu is a homogeneous function of X12 of degree 2. It implies 
that 

(rf - 2 + X12 ■ Vu)H,2 = id + 2)Hu + [ . . . ]i2 • 

Comparing this with the relation (|32D, we obtain from eq. (p9| ) the value of the coefficient 
a: 



id-l){d + 2) 
Next, it follows from the relation (^) that 



[ • • • J12^ 
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Comparison of eqs. p8| ) and ( |3lD gives the value of the coefficient b: 



1 / N-2 d + N-2' 

"d V (d- l)(rf+2) ^ 2 



This completely determines H12 up to terms 



12- 



Finally, in order to find the anomalous dimension F^q, we recall that the matrix of anoma- 
lous dimensions is found by looking at J2jytkHjk, cf. eq. (p7[). F°q is obtained by projecting 
the relation ( pTj) on Eq using the triangular structure of the transformation F. Gathering all 
the terms in eq. (p3|), we obtain after a simple algebra 



^ ^'^ " " 2(d + 2) ^ 1^ 

j^k ^ya-r ^) pairings {{;_,;+)} (/_^;^) 



l<i_<i+<iV 

iV(rf + iV) 



2(d + 2) 



^0 + [ 



Thus 



^ Ar(rf + N) ^_N_ NjN - 2) 

° 2(d + 2) 2 2(rf+2) ^ ' 



which via eq. (pSf) leads to the claimed value (Bl) of the anomalous exponent pn- 



5 Dissipation field 

Let us sketch the argument for the equality of two definitions (0J^) of the dissipation field 
e{x) at = 0. First, we shall retrace the self-consistency arguments about the short distance 
behavior of the correlation functions [Q, ^. These go as follows. By separating terms, the 
Hopf identity ( |T0|) may be rewritten in the form: 

(^-z/Ai - + \D{d - l)A^(i,2)) (T(a;i)T(x2)T(x3) . • • T{xn)) 
= E [^^3 - \D{d - l)(A^(i,,) + A^(2,,)) (T(xi) . . . T{x^)) 

j=3 ^ ^ 

-'-Did-1) J2 -Mo»(T(xi)...T(a;jv)) + J2Cix,k/L) {Tix,). .T{x^)) (38) 

3<j<k<N j<k 3 k 

where the 2-point operator ^D{d — l)Ai(^j^k) = d°'^{xjk)dx"d^i3. In variables x = and 

^ k 

y = X12 the l.h.s. of eq. (^) becomes 

(^-2uAy - dycdy, - iz/A,. + 9^.9,,) {T{x,) . . . T{xm)) . 

Eq. (|38D , with the use of the latter decomposition and of the relation d^a = d^^ + d^^^ = 

N 

— J2 dx'^ 1 allows to write 

(-2Z/A, - dy^dy,) {T{x + \y)T{x - \y)T{x^) . . . T{xn)) = R (39) 
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where i? is a combination of terms involving only Xj-derivatives and at most first y-derivatives 
of (T(a; + iy)T{x — ^y)T{xz) ■ ■ .T{xn))- Let us assume that the limit when Xi2 — >■ of 
(^{xi)T{x2)T{x^) . . . T{xn)) and of {(VT){xi)T{x2)T{x3) . . . T{xn)) and of their derivatives 
over Xj, j > 3, exists uniformly in small u > (for separated x, X3, . . .xn)- Then, as in the 
analysis of the 2-point function (with anisotropic forcing), one infers from the eq. ( |39D that 

{T{x + \y)T{x - \y)T{x:,) . . . T(x^)) = cM\^ + \D{d - l)\y\^)-^ 

+ zero modes of {—2uAy — d°''^{y) dyadyp) + error (40) 

with the coefficients depending on x,X3, . . .xn and the error more regular when u —>■ and 
y ^ 0. The zero modes contain a polynomial of the first order in y. Of the remaining zero 
modes the most dangerous one comes from the angular momentum 2 sector and it behaves 
like 0{\y\'^^) for u < hD{d-l)\y\^, where ^2 = i(-rf + 2 - ^ + y^(rf - 2 + 0^ + 8d) > 2-^, 
and like 0(|?/p) for ^D{d — l)\y\^ ^ i^- All such terms and their first y-derivatives have 
limits when y — > uniformly in small u. As we see, our assumptions about the correlators of 
T are at least self-consistent. They are confirmed by our 0{C,) computation of the structure 
functions. Indeed, at = and for large L the structure functions receive the dominant 
contribution from the zero modes Fq of Ai^ which behave like ^O^yl"^) log \y \ modulo a first 
order polynomial in y and 0(^^) terms, in agreement with the above analysis. Note that the 
^0(|?/P) log \y\ contribution to Fq is not rotationally invariant in y: it receives contributions 
from both the (9(|?/p~'^) and the 0{\y\°'^) angular momentum 2 terms in {T{xi) . . .T{xn))- 

Let us use our self-consistent assumptions about {T{xi) . . .T{x2)) in a version of eq. 

(2uVi-V2 + d"^{xu)d,<^d^,){T{x,)T{x2)T{x;)...T{xN)) = R' . (41) 

Expression R' involves only terms with at most one derivative over xi or X2- Therefore 
lim lim R' should exist and be equal to lim lim R'. The same limits applied to the l.h.s. 

of (^ID give, depending on the order, the definitions (0) or (H) of the dissipation field insertion 
e{x), provided that ^ > 0. Indeed, under first xu — ^ and then u ^ limits the (i(xi2)ViV2 
term disappears due to vanishing of d{x) at zero while sending u —* before the xu — > 
limit kills the z/Vi ■ V2 contribution. Hence the equivalence of two definitions for < ^ < 1. 

By similar arguments, all three limits lim, lim and lim commute in the action on R'. 
Applying them on the l.h.s. of eq. (PI), we infer that at z/ = 

lim e{x) = '-Did -l)iVTix) f (42) 

5— »o ^ 

and it describes the dissipation field of the scalar T diffusing with the diffusion constant 
^D{d — 1) and dissipating energy on long scales. Note that the r.h.s. may be viewed as 
a direct application of the second definition (§) at ^ = whereas the application of the 
first one (|^) would give e{x) = 0: the equivalence of the definitions breaks down at ^ = 0. 
At z/ = 0, the x' —>■ X and ^ — > limits do not commute for {'VT){x'){VT){x) although 
they do commute for T{x')T{x) or for {'VT){x')T{x). This is due to the disappearance 
of the distinction between the dissipative and the inertial range behavior at ,^ = 0. A 
straightforward calculation employing the relation (^) shows that at non-coinciding points 

lim lim {e{xi),...,e{xn)y = 2"-i(n - 1)! rf^'^e" . (43) 
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In particular, field e{x) becomes constant in space at ^ = and L = cxd, in agreement with 
the physical picture of dissipation becoming a large scale phenomenon when ^ — > 0. 

The inertial range decay (^) follows from eq. ( |T6|) with the use of the definition (Q) of the 
dissipation field and of the fact that (e(x))^ = gives for large L a subdominant contribution 
to (e(x), e(0))^. From eq. ( ^31) we infer that the proportionality constant in (H) is equal to 
^ + 0{^). Similarly, the mixed correlation functions (e(xi) . . . e(x„)T(?/i) . . .T{ym)) scale 
with the infrared cutoff as Lp^"+"^ and with positions with exponent (2 — — p2n+m- 



6 Equations for structure functions 

Much of the past attempts to understand the behavior of the structure functions Sn was 
based on the differential equations satisfied by them These equations, may be ob- 
tained from the iV-point function equation ([lOD in the following way. Let Sj{x,y) denote 
the difference operator acting on functions of N variables /(xi, . . . , xn) by subtracting their 
values at Xj = x and Xj = y. Sj{x,y) commute for different j and {(T{x) — T{y))^) = 
Y[j Sj{x, y) {T{xi) . . . T{x]\f)) . Apphcation of Yij y) to eq. ([10|) results in the identity: 

- d''P{x)d,.d,0 Sn{x) + N{N -1){C{1) - C{Q)) Sn-2{x) = Jn{x) (44) 
with the dissipative contribution 

Mx) = 2z/iV((AT)(x)(r(x)-r(0))M- (45) 

Alternatively, eq. (0) may be obtained directly from the basic stochastic differential equation 
(IID. Despite prefactor z/, the term Jjv does not vanish when z/ — due to the dissipative 
anomaly. Indeed, coefficient u may be absorbed into the insertions of the dissipation field: 

Jjv(x) = 2uASn{x) - 2N{N -l){e{x){T{x) -T{0)f-^) 

— > -2Ar(Ar-l)(e(a:)(T(a:)-T(0))Ml.=o. (46) 

The mutual non-sing ularity of e(x) and {{T{x') - T{0))^ at X = x' and u = 0, which has 
been assumed in the last expression for Jjy, may be checked by a self-consistent analysis or 
directly for the perturbative solution. 

Our result (|) about the asymptotics of the structure functions and eq. (|4^ ) imply that 
at u = and for large L 

with bN = D{d-1) ((2 - Of - Pn) [d+{2- 0^ - Pn) aN ■ Note that the last expres- 
sion may be rewritten as 

(48) 

This relation may be confirmed by a direct calculation of Jjv to order 0{^) from the dominant 
zero mode Fq contribution to the correlation functions. 
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In the inspiring paper 0, Kraichnan attempted to obtain anomalous exponents from 
eq. (11) by assuming a relation similar to (||) but with = ((^-O f -pn ^ -p«) 



replaced by y. His assumption led upon insertion to (|^) to the quadratic equation for the 
scaling dimension C,]si = {2 — — 



CN{CN + d-2 + i) = {2-i)d 



N 



whose solution gave the anomalous exponents p^v- Note that the replacement of the factor 

a2bN 
ajvb2 



N 



on the r.h.s. of Kraichnan's equation for by leads instead to a tautological identity. 



Kraichnan's values of p^v, unlike the ones obtained in the present work, do not vanish 
at ^ = 0. The latter might seem strange in view of the fact that T becomes a Gaussian 
field at ^ = and the Ansatz of was fit with the Gaussian calculation. The latter was 
based, however, on the definition ( ^5]) applied directly in the Gaussian case whereas at ^ = 
one should use for Jn the expression (|46|) with e{x) given by the equation (|^). The latter 
calculation agrees, of course, with eq. (pH]): 

Jn{x)\^=, = -b%\xf-' (49) 

where 6^ = D{d — l)N{d + N ~ 2)a^, see eq. (0). In this limit, the differential equation 
reduces to —D{d — 1) ASn\^=o = Jn\^=o- 



Above, we have studied the inertial range behavior of the structure functions. More gen- 
eral objects to study are the joint probability distribution functions (p.d.f.'s) P/v(Ti, . . . , T/v; 
Xi, . . . ,xn) of the scalar whose moments give the equal time correlation functions. In par- 
ticular, the structure functions (Oj iT{xj) — T{yj))) are special moments of the p.d.f.'s 

QNiTi, ...,Tn;xi,...,Xn) = J PNiTi + r, . . . , T^v + r ; Xi, . . . , Xat) dr (50) 
which are translational invariant in the T- variables. The generating function of S'at's, 

Z{\-x) = (e*^(^(-)-no))^ ^ f e^^^Q{T;x)dT (51) 



is a Fourier transform of the p.d.f. Q(Ti — T2; 0:12) = Q2iTi,T2; xi, X2), 

In a recent paper on the Burgers equation, Polyakov has argued that the structure 
function p.d.f.'s exhibit a universal inertial range behavior for (translating his statements 
to the passive scalar case) Tj <^ T^ms where Tj-ms = i/(T(0)^) = (9(L^~^/^), see eq. (|14[). 



Polyakov's analysis was based on postulating an operator product expansion which allows 
to close the resummed version of eq. (H4|) 



- ci"^(x) d^cd^p Z{\; x) + ( C(0) - C( J) ) Z(A; x) = J(A; x) (52) 

by expressing its r.h.s. J(A; x) = J2n ■^^^rp'^Afl^^) = 2A^(e(x) e*'^*^'^^^)~'^^°^)) again in terms of 
Z{X;x). The resulting equation for Z may be reduced to an ordinary differential equation 
by imposing the scaling relation 

(2x- V- (2-0A9a) ^(A;x) = 0, (53) 
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i.e. by postulating that Z{\; x) depends on A^|xp"^. Note that a strictly scaling solution for 
Z(A; x) implies either the Kolmogorov scaling or divergence of the structure functions. For 
example, the resummation of the expressions ( p^Tf ) for the ^ = 0, z/ = 0, L = oo structure 
functions gives the Gaussian generating function 



Z{\;x)\^=o = exp 



Dd{d-l 



\2| |2 

■A \x\ 



(54) 



which scales, in accordance with the normal scaling of the ^ = structure functions. It 
corresponds to the Gaussian p.d.f. 



Q{T-x) 



'Dd{d-l] 



1/2 



— exp 

\x\ 



Dd{d-l) 



4e 



(55) 



A straightforward check shows that at ,^ = the dissipative term J may be simply expressed 
in terms of Z itself: J = 2lX^{l + i Xdx) Z. 

Our small ^ analysis of the passive scalar does not allow us to confirm or to infirm 
Polyakov's picture since the individual structure functions that we study probe the small A 
behavior of Z i.e. the large T behavior of Q{T). What follows from it is that the large T 
tails of the p.d.f. 's violate scaling for ^ > 0. In particular, the solution (|],^ for Sn^s leads 
to the relation 



2x- V - (2 



d + 2' d + 2 

and to the asymptotic scale L dependence 



{\dxf]Z{\;x) 



(56) 



Ldr- 



Xdx{Xd,-2)\Z{X;x) = 0{e) 



2{d + 2) 

consistent, of course, with the exact overall scaling property of the problem: 

{2x ■ V + 2 LOl - {2 ~ ^dx) Z{X; x) = 0. 



(57) 



It would be interesting to recover the tail of the structure-function p.d.f. Q in the order 
0{^). This would require the knowledge of the 0{C,) contributions to the non-universal 
amplitudes in the relation (|^) which we have not computed. It is possible that they may 
be found by a perturbative analysis of instanton contributions to the functional integral for 
the dynamical scalar correlators ll^]. We expect that the refined perturbative analysis in ^ 
will allow a better control of both the structure functions and the corresponding p.d.f. 's. 
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